Condensation of Anyons in Frustrated Quantum Magnets 
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We derive the exact ground space of a family of spin- 1/2 Heisenberg chains with uniaxial exchange anisotropy 
(XXZ) and interactions between nearest and next-nearest-neighbor spins. The Hamiltonian family, H c ft(Q), is 
parametrized by a single variable Q. By using a generalized Jordan-Wigner transformation that maps spins into 
anyons, we show that the exact ground states of H C «{Q) correspond to a condensation of anyons with statistical 
phase 4> — —4Q. We also provide matrix-product state representations of some ground states that allow for the 
efficient computation of spin-spin correlation functions. 

PACS numbers: 



One-dimensional quantum magnets can realize exotic states 
of matter such as Luttinger liquids fl] |2], valence bond 
solids J3JIH, and spin supersolids [5|. A unique feature of 
these systems is that transmutations of particle statistics pre- 
serve the range of the interactions. This is the main reason 
behind the success of spin-fermion transformations for study- 
ing one-dimensional spin models [6-8|. A simple general- 
ization of the Jordan-Wigner transformation [ 6 1 allows us to 
map spins into anyons (particles that generalize the concept 
of fermions and bosons). By exploiting this generalization, 
we derive exact ground states of spin- 1/2 frustrated magnets 
and provide an efficient method for computing spin-spin cor- 
relation functions. Remarkably, the resulting ground states 
are anyon condensates that spontaneously break the particle- 
number conservation symmetry. In contrast to the familiar 
Bose-Einstein condensates (BECs), the condensed particles 
satisfy anyonic statistics. 

BECs are ubiquitous in many-body physics. Atomic gases, 
certain quantum magnets, and superconductors can be de- 
scribed as interacting gases of bosons that become BECs at 
low enough temperatures l9l- [T3l . A BEC is characterized by 
the spontaneous breaking of the symmetry associated with the 
boson-number conservation. This symmetry breaking comes 
along with a macroscopic occupation of the single-particle 
ground state, suggesting that the notion of BEC cannot be gen- 
eralized to other particle statistics. Indeed, anyons, which are 
identical particles obeying arbitrary statistics [14], are sub- 
ject to a generalized Pauli exclusion principle that limits the 
occupancy of any given state. Therefore, the only way of gen- 
eralizing the notion of BECs to anyon condensates (ACs) is 
by taking an adequate limit over the particle statistics that sur- 
vives the exclusion principle. Like BECs, ACs spontaneously 
break the particle-number conservation symmetry. 

Just as BECs can only exist at zero-temperature in low- 
dimensional systems with short range interactions Q3], we 
will demonstrate that ACs are the ground states of a family 
of one-dimensional frustrated quantum magnets. The anyonic 
statistics is fixed by a phase <j> that their wave function acquires 
after exchanging two anyons: vf^ri,^) = e l ^ 1 ^(r2, ri). This 
condition implies that anyons are in general "non-local" par- 
ticles that interpolate continuously between bosons (<f> = 0) 
and fermions (<f> = tt). 



We will start by considering the dimerized S — 1/2 spin 
system of Fig. [TJ A. Pairs of spin-1/2 particles in the same 
dimer interact via dominant antiferromagnetic exchange (dou- 
ble lines). Weaker inter-dimer exchange interactions are rep- 
resented by single lines. In presence of an external field h, the 
low-energy single-dimer states are the triplet |ff) (s z = 1/2) 
and the singlet |fJ,) - \H) O z = -1/2). The low-energy 
effective spin-1/2 XXZ model, 7i c s , that results from project- 
ing the original Hamiltonian into the low-energy subspace is 
EMU (see Fig. [T): 
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The eigenstates of s| represent the singlet and the triplet of 
the j-th dimer Hamiltonian. 

In this Letter we derive the exact ground states of H c s(Q) 
for A„ = cos(^Q) and J\ = — 4J 2 cosQ, with < 
Q < tt. HQ = 7r/2 the model reduces to two decoupled 
ferromagnetic Heisenberg chains whose exact solutions are 
known [20|. Consequently, if the number of effective spins 
in the ladder is L, the ground-space dimension is (L/2 + l) 2 
or [L + 3)(L + l)/4 for even or odd L. We will prove that 
the ground-space dimension is independent of Q for a special 
choice of open boundary conditions. 

Some of the ground states are product states \sg(Q)), in 
which the spins are aligned according to a canted spiral of 
wave vector Q and arbitrary canting angle 8 [21 J (see Fig. [2] 
A). \sg(—Q)) are also ground states due to spatial inversion 
symmetry: H s s(Q) — H e s(—Q)- The long-range spiral or- 
der is manifested by the spatial dependence of spin-spin corre- 
lations that oscillate with constant amplitude and wave-vector 
±Q(seeFig.|2]B). 

States with spiral order can be represented as BECs by 
mapping spin-1/2 operators into hard-core bosons via the 



st. The 



Matsubara-Matsuda (MM) transformation (221 : b 
operator creates a hard-core boson at the j-th site of the lad- 
der. The identity (frj) 2 =0 excludes states with more than one 
boson per site. Since H e s(Q) is translationally invariant, it is 
useful to introduce the operators b\ = Ylj=i exp(ifcj)&| / \[L 
that create bosons of momentum k. Contrary to real space 
bosons, bp occupancy in momentum space is not subject to 




S = 1/2 



FIG. 1: Frustrated quantum magnet. (A) Original spin- 1/2 system: 
pairs of spin- 1/2 particles in the same dimer interact via dominant 
antiferromagnetic exchange (double lines). Weaker inter-dimer ex- 
change interactions are represented by single lines. (B) Effective 
spin- 1/2 XXZ model, H c a, that results from projecting the original 
Hamiltonian into the low-energy subspace. 
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FIG. 2: Canted spiral order. (A) Spin- 1/2 representation of the spiral 
states \sg (Q)) for Q — 2n/5. (B) Spin-spin correlation as a function 
of distance r. 

exclusion. The BEC states 

\M±Q)) cx (b{ Q ) n \0) (1) 

span the ground space with spiral order [21]. |0) is the vac- 
uum or empty state that corresponds to the fully-polarized 
state \l ... I) in the spin language and n is the number of 
single bosons with momentum ±Q. \ip n (Q)) and \ip n {— Q)) 
are usually non-orthogonal but linearly independent. Finding 
all BEC states for Q = requires to consider adequate lin- 
ear combinations of \ip n {Q)) an d iV'nC - Q)) before taking the 
limit Q — > 0. The single-boson excited eigenstates, b\.\0), 
have energy ujk — 2 J 2 (cos k — cos Q) 2 [21 1. 

The MM transformation provides then a simple alterna- 
tive representation of the spiral ground states in terms of 
hard-core bosons. Remarkably, a generalization of the par- 
ticle statistics will allow us to find exact expressions for all 
ground states of H c s(Q)- The generalized JW transformation 



a}- = exp{i0^ i<:) (sf + l/2)}s^ maps spins into anyons, 

where the operator a - creates a hard-core anyon [(a]) 2 = 0] 
at the j-th site of the ladder. The MM and JW transformations 
are recovered for <fi = (hard-core bosons) and <fi = it (spin- 
less fermions), respectively. The anyonic statistics becomes 
manifest in the commutation relations a)a\ = e-^aja], 
when j < I. The states 

|VW*(Q))= ton c (4)"(al Q r|0) , (2) 

with < n,m < L/2, provide a natural generalization of 
the BEC states in Eq. ([T}. a/ k = ^ J=1 exp(zfc i j)aj/v / i is 
the operator that creates an anyon with momentum k and 
is the normalization constant. The limit in Eq. (|2| is neces- 
sary to avoid a trivial cancellation due to a generalized Pauli 
exclusion principle: (ajj p = for 4> = 2ir/p and p an in- 
teger greater than 1. The statistical phase does not play a 
role if n — or m — because lirri0^_4Q c^ag)™! ) (yi 
(6g) ra |0). However, the same limiting procedure leads to 
states that are qualitatively different from a double-Q BEC, 
(6^)"(&^ Q ) m |0), when n,m ^ 0. 

\ipn,m{Q)) span the full ground space of H e s(Q) - see be- 
low. As n + m is not preserved, the particle-number symmetry 
is spontaneously broken. Therefore, we refer to the physical 
phenomenon characterized by such states as anyonic conden- 
sation and to \ipn,m(Q)) as AC states. The physical properties 
of ACs are better understood if we first consider the action of 

. Besides polarizing the j-th spin, such operator creates a 
"phase kink": it rotates the spins on the left of j by an an- 
gle <j) along the z axis. The phase kink gives rise to the any- 
onic statistics that is absent in the MM transformation. Conse- 
quently, the ground state ag|sg(— Q)) corresponds to a phase 
soliton attached to a boson that propagate with momentum Q 
across the spiral - see Fig. [3] A. In general, the ground states 
(ag) ra |se(— Q)) will correspond to a condensation of solitons. 

To demonstrate that \ipn,m{Q)) is a ground state of 
M e fi(Q)> we rewrite the H c s{Q) as 

L-l 

n cS (Q) = £(j 2 + j 2 /(8J 2 ))n,(Q) , (3) 

3=2 

where the n^Q) are projectors acting on triangular plaque- 
ttes ETI . These projectors are 

n J -(Q) = i4xeJi + i4)^i ! 

where oc (2 cos Qs+_ r s+ +1 -s+s+ +1 -s+_ 1 ,s+)||||) and 

oc (2cosQs+ - s+_ x - s1 +1 )|44-i) are plaquette states. 
Note that Eq. ([3]) refers to a particular choice of open boundary 
conditions. Since the fully-polarized spin states |J, . . . I) an d 
|t ■■■ t) are orthogonal to |£J), they are zero-energy ground 
states of Wctf(Q)- This implies that H c s(Q) is a so-called 
parent Hamiltonian and any other ground state must have zero 
energy. 



2 



The AC states satisfy 
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jl < • • • < 3: 
k < ... < k 
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The limit in Eq. |4} removes the prefactors that arise from ex- 
changing two anyons with same momentum. This is necessary 
to avoid a trivial cancellation of (oj_q)™ caused by the gener- 
alized exclusion principle. In particular, the limit implies that 

|^n,o(Q)) oc (¥ Q ) n \0) and |^ ,m(Q)) oc (P_ Q ) m \0). Let 
v p = {p — l,p,p + 1}. Then, if j q ^ v p and l r ^ v p for all 



n P (g) n a ] q n a |j0) = n «j, n ^wi") = °- 
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We now consider the terms in Eq. Q for which a single ji £ 
v p and /,. ^ v p for all r. We fix all the indexes except for j, 
and consider the sum over the terms that have ji in v p . The 
resulting state satisfies 

p+l n m 

n P (Q) E U eiQj * a lU e ~ iQlra V ) * 

ji—p—1 q—1 r—1 

n ^'"i, n e_4Qira t n p(Q) e eiQii 4 1 ) = ° ■ 



j-=l 



This equality follows from Ejjt=p-i e'^'ajj . |0) being or- 
thogonal to A similar reasoning applies if a single l r 
is in Up and j q ^ i!p for all q. Finally, we analyze the case 
for which only two indexes of {j q } q U {l r } r are in v p . If the 
two indexes are either in {j q } q or {/ r } r , we can use a similar 
argument as the one above to show that H p (Q) acts trivially 
in the particular sum of states. We then analyze the case for 
which one of the indexes in v p is in {j q } q and the other one is 
in {l r } r . We label these indexes by ji and 1^, respectively. If 
we fix all the remaining indexes and consider the sum of terms 
over the six possible values of (ji, If,,), the resulting state sat- 
isfies 

n tci 

n P (Q) e n^Hn^Hi )* 
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The state £V i h ) e lQ( - ji ' lh) a^a\ h \0) is orthogonal to 
when 

e iQ + ei (3Q+0) _ cos q (1 + e ,(4Q+0) ) = ) 

where <f> is the anyonic statistical phase. This condition is sat- 
isfied for (f> — — AQ. 



For even L, T-L e g(ir/2) corresponds to two decoupled ferro- 
magnetic Heisenberg chains of length L/2. The ground space 
is spanned by the linearly independent states \^ n<m (ix/2)), 
with < n,m < L/2, and the ground space dimension is 
(L/2 + l) 2 . Because the AC states are analytic in Q, the 
set {\ip n ,m(Q)}n,m remains linearly independent and the di- 
mension of the corresponding subspace does not change with 
Q. (Numerical calculations suggest that these are all possible 
ground states.) The states \ip n ,m(Q)) an d IV , m,n( — Q)) be- 
come equal for Q — > 0. However, it is possible to obtain two 
linearly independent states by taking a proper limit of linear 
combinations of \ip n ,m{Q) an d \^m,n{Q)- For instance, 



lim 

Q-yO 



E Ul + ---+3n)b] 1 ...bl\i 
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J'l<-<J'» 



is also a ground state of H c s(0). 




FIG. 3: Phase soliton. (A) A topological defect results from the 
creation of an anyon on the spiral state, i.e. aAs$(—Q)) (Fig. 2 
A). (B) Exponential decay of spin-spin correlations as a function oi 
distance r for an AC. This correlation function was obtained by using 
a matrix-product state representation. The Fourier transform (FT) of 
the correlation function is peaked at the momentum of the AC. c > 0, 
Q = 27r/5, andL = 100. 



The difference between ACs and BECs is evident from 
their single-particle correlation functions. The two-point 
correlator exhibits long-range order, lim,._ i . 00 \ (sjSj +r )\ = 

linWoo K&j-&7 +r }| ^ 0, for a BEC (Fig.[l]B). In the anyonic 

maps to exp{i^JVj r }ata. +r , 

counts the number of anyons 



language, the operator SjSj +r 



where iV--, 



, / i+r a \ a l 



between sites j and j + r. Then, the large fluctuations of Nj r 
present in an AC translate into large fluctuations of the rela- 
tive azimutal phase between spins j and j + r— see Fig. [3] A. 
If we assume that the probabilities of having anyons in dif- 
ferent sites are roughly independent of each other, we obtain 
(s+sj +r ) ~ e- iQr (l - p(l - e^)) r for large r and p = n/L. 
This simple analysis suggests that the two-point correlator in 
the ACs decays exponentially in r - see Fig.[3]B. 
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Spin-spin correlations can be efficiently computed by in- 
troducing matrix-product state representations for the ground 
states that are invariant under translation [23 24]. These states 
are represented as 

MQ)) ex Y, tv[A^(Q)...A^(Q)]\a 1 ...a L ),(5) 



and u) = cxp(i27r/p). Few-spin correlators can be numeri- 
cally obtained with a low computational cost that is polyno- 
mial in p 11241 . For instance, 

WQ)\4-sj+MQ)) = tr x (AHQ) ® (^(Q)Dx 
xE r ~ l x ® (AHQ)T) x £ L - r ] , (9) 



where each A"' (Q) is a matrix associated with the j-th spin 
in the ladder and o,- e i}. Correlation functions can be 
computed analytically or numerically in this representation 
by a contraction of a tensor network 11251 . The translational- 
invariant ground states of the matrix-product form of Eq. <|3j 
satisfy 



(6) 



The action of U p (Q) on \ip(Q)) outputs another matrix- 
product state whose three-spin tensor, associated with v p , is 

B t p-i< t p> t p+ 1 (Q) = (7) 
H A ap - 1 (Q)A a "{Q)A a ^ 1 (Q)x 

x (t p -i,t p , r p +i|Hp(Q)|op_x, Op, ct p +i), 

and r, o € {til}- Equation |7]i follows from the contrac- 
tion of a tensor network GSJ. B Tp_1 ' Tp,Tp + 1 (Q) = for all 
T p-ii r p> T p+i is a sufficient condition to satisfy Eq. (|6j. It 
follows from the definition of ) that: 

^(Q)^(Q)^(Q) + ^(Q)^(Q)^(Q)- 
- 2cos(Q)^(Q)At(Q)^(Q) - , 

A±{Q)A\Q)tf{Q) + A^(Q)A^Q)AHQ)- 

-2cos(Q)A^(Q)AHQ)A^(Q)=0 . (8) 

These equations are satisfied if A^(Q)A^(Q) = 
e ±lQ A^(Q)A^(Q). Equations ^ are also satisfied for 
Q = 2ir/3 if A^'^(2ir/3) are the matrices that correspond to a 
ground state of the well-known Majumdar-Ghosh model 11261 . 

Solutions for Q = 2ir/p, with integer p, are given by 
A^(Q) = Z p and A±(Q) = (aX p + bX\)Z\. The com- 
plex constants a and b are arbitrary and the corresponding 
states span the full translational-invariant ground subspace. 
Z p and X p are p-dimensional unitary matrices that satisfy 
ZpXp — e~ % QXpZ p . This condition resembles the anyonic 
nature of our solutions whose statistical phase is <fi = — 4Q. 
For example, Z p is the diagonal matrix and X p is the cyclic 
permutation that generate the finite Heisenberg or generalized 
Pauli group of order p 2 



Z P = 



(1 ••• \ 
u ••• 



\0 



/0 1 ••• 0\ 
1 ••• 



) Xp — 



\1 



and E = A*(Q) ® (A^(Q))* + A±{Q) ® (A+(Q))*. The 
normalization constant is (■0((5)|'(/'( ( 3)) = tr[25 i ]. We used 
Eq. |9} and Gram-Schmidt orthogonalization to compute the 
exact spin-spin correlator in Fig.[3]B and confirm the expected 
exponential decay for ACs. 

Finally, it is interesting to note that the particular ground 
state solutions ([TJ imply that the scattering length of two 
bosons with momenta k and k' goes asymptotically to zero for 
k,k' —> Q or k, k' —> —Q. However, the boson-boson scatter- 
ing length remains finite if k — >• Q and k! — > — Q. In contrast, 
our general ground state solutions given in Eq. (|2) imply that 
the scattering length of two anyons with momenta k and k' 
goes asymptotically to zero in both cases: (k, k') — > (Q, Q) 
or (k,k') -> (~Q,-Q) and(k,k') ->■ (Q, -Q). This im- 
plies that the low-energy excitations on top of the condensate 
are quasi-free anyons. 

In summary, we have introduced the notion of anyon con- 
densates as a generalization of BECs. By using a natural ex- 
tension of the Jordan- Wigner transformation, we have shown 
that that the AC states (|2]i are the exact ground states of a 
family one-dimensional XXZ magnets with nearest and next- 
nearest-neighbor exchange interactions. ACs spontaneously 
break the particle-number conservation symmetry and corre- 
spond to a condensation of topological defects. We have also 
provided a matrix-product representation for the translational- 
invariant ground states, which allows for the efficient compu- 
tation of correlation functions. These states also satisfy an 
area law, which is rather common in frustrated systems, and 
results in an upper bound on the entanglement entropy 11281 . 
The exponential decay obtained for the single-particle correla- 
tor in the AC states is a natural consequence of the topological 
nature of anyons. Since each anyon is a bound state of a boson 
and a phase soliton, the large fluctuations in the anyon number 
that characterize an AC state lead to large phase fluctuations. 
This behavior is qualitatively different from the bosonic limit 
(BEC) in which particle number and phase are conjugate vari- 
ables and large number fluctuations lead to long-range phase 
coherence. 

Work at LANL was performed under the auspices of the 
U.S. DOE contract No. DE-AC52-06NA25396 through the 
LDRD program. 
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